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Abstract 

In this note we concern with the wave maps from the Lorentzian manifold with the periodic in 
time metric into the Riemannian manifold, which belongs to the one-parameter family of Riemannian 
manifolds. That family contains as a special case the Poincare upper half-plane model. Our interest 
to such maps is motivated with some particular type of the Robertson- Walker spacetime arising in 
the cosmology. We show that small periodic in time perturbation of the Minkowski metric generates 
parametric resonance phenomenon. We prove that, the global in time solvability in the neighborhood of 
constant solutions is not a stable property of the wave maps. 
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1 Introduction 

In this note we prove instability, under periodic in time perturbation of the Minkowski metric, of the global 
in time solvability of the wave map equation in the neighborhood of stationary solutions. 

Let (V, g) be Lorentzian and (M, h) be Riemannian manifolds of dimensions n + 1 and d, respectively. 
Let u be a continuous mapping from V into M: 

u : V — > M. 

The wave map is the mapping u that satisfy the Euler-Lagrange equations for the Lagrangian 

C{u) '•= J ^(du,du) T , V(S)u ~.i TM dV 

with dV = y/gdxi A ... A dx n +\ in local coordinates, being the volume form of V. From now on the Einstein's 
summation convention is in force. Written in the local coordinates on V and M the Euler-Lagrange equations 
read 

g afi (d 2 afi u A - T x a pd x u A + Ti c d a u B d p u c ) = . (1) 

Here and Tg C denote the components of the Riemannian connections of g and h, respectively. The 
wave map equation is invariant under isometries of (V, g) and (M, h). 

We are interested in the Cauchy problem associated with ([1]). More precisely, for given initial data 
(u(0), dtu(0)) : So — > M x TM at time t = 0, we look for a wave map u extending these globally in time. 
First, consider the case of (V, g) being the standard Minkowski space equipped with metric 



g°P = diag(-l,l, !,...,!) • 



(2) 



In connection with this metric we recall conjecture of Klainerman and related known results for the Poincare 
upper half-plane model (the standard hyperbolic plane). 

Conjecture (Klainerman |12j ) Let (M 2 ,h) be the standard hyperbolic plane. Then classical wave maps 
originating on R 2+1 exist for arbitrary smooth initial data. 

The following partial result towards the conjecture has been established in [12]. 

Theorem 1.1 (Kreiger [12]) Let (M 2 ,h), H 2 := {(u 1 ,^) € R 2 | u 2 > 0},, be the standard hyperbolic plane 
with the metric tensor hndu l du-' — , „ Nn ((du 1 ) 2 + (du 2 ) 2 ). Then given initial data u[0] : {0} x R 2 — > 



(u 2 ) 2 

x TH 2 which are sufficiently small in the sense that 



E 



l 2 



d a u 2 



for suitably small e > 0, there exists a classical wave map from 



dx < e (3) 
1 to H 2 extending these globally in time. 



In particular, the global wave map exists for any small, in the sense of ([3|), compactly supported perturbations 
of the initial data of the constant wave map, which has vanishing integral ([3]). In other words, the global 
solvability of the wave map equation is a stable property in the small neighborhood of the global constant 
solutions. 

The answer to the Klainerman's conjecture as well as the scattering result for the wave map are given 
in |13j . In particular, it is proved in |13j that if M is a hyperbolic Riemann surface, and initial data 
(u(0), dtu(0)) : So — > M x TM are smooth and u(0) = const, dtu(0) — outside of some compact set, 
then the wave map evolution u of these data as a map ]R 2+1 — !• M exists globally as a smooth function. 

In this paper we are interested in the case of the Riemannian manifold (M, h) which belongs to one- 
parameter family of manifolds containing the Euclidean half-space and the Poincare upper half-plane model 
(M 2 ,h). In fact, that family consists of the Riemannian manifolds, which are the half-plane {(u 1 ,u 2 ) £ 



> 0} equipped with the metric hijdu l du 3 



1 



(u 2 ) 1 



[(du 1 ) 2 + (du 2 ) 2 ), where the parameter / is a real 



number. For I = the metric is Euclidean, while for I — 2 it is the metric of the standard hyperbolic plane. 
Those are the only two manifolds of this family which have constant curvature. 

In the present paper we examine the stability of the global solvability of the wave map equation with 
respect to the perturbation of the metric g. First, we prove that the only stationary solutions of the 
equation (TTJ) are the constant solutions. Then, we show that the global in time solvability can be destroyed 
by parametric resonance phenomena. (For the scalar quasilinear wave equation it was proved in |23].|24j.) 
Thus, the small data global solvability is unstable with respect to the arbitrary small periodic perturbation of 
the metric tensor g. More precisely, we prove that the local solution obtained by small, for given G S N in 
the sense of the following integral 




d^u 1 



(u 2 ) 1 ' 2 



d"<i 



(u 2 ) 1 / 2 



dx, 



(4) 



smooth compactly supported perturbations of the initial data of the global solution, in general, cannot be 
extended globally in time. For the parametric resonance phenomena in the scalar wave map-type hyperbolic 
equations see [25] and references therein. Then, according to [32] (see also references therein) the parametric 
resonance phenomena in the linear scalar wave equations can be localized in the space. 

The Cauchy problem for the wave maps in the perturbed Minkowski spacetime is considered in i3: . More 
precisely, assume that V = S x R, with S an n-dimensional orientable smooth manifold, and let g be a 
Robertson- Walker metric g = —dt 2 + R 2 (t)a, where a = o~ij dx 1 dx 3 is given smooth time independent metric 
on S, with non-zero injectivity radius. The Christoffel symbols of the Robertson- Walker metric are 



1 00 — 1 j0 ^ 1 00 — u i 



Fq^- = R ^R'Sj, T®j — RR'o~ij, R' : — 



clR 



]li l]ln - uu -ju - uu -i - uj j ' - 13 v -- ■ ^ 

where 7 denotes the Christoffel symbols in the metric a. The following result is known. 
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Theorem 1.2 (Choquct-Bruhat [3]) Let (S xM, g) be a Robertson-Walker expanding universe with the metric 



<) 



dt 2 + R 2 (t)a, 



(5) 



with (S*, er) a smooth Reimanian manifold with non-zero injectivity radius, of dimension n < 3, and R a 
positive increasing function oft such that 1/R(t) is integrable on \to,oo). 

Let (M,h) be a proper Riemannian manifold regularly embedded in M. N such that Riem(h) is uniformly 
bounded. 

Then there exists a global wave map from (S x [to,co),g) into (M,h) taking Cauchy data if, ip with Dip 
and ip in Hi if the integral ofl/R(t) on [to,oo) is less than some corresponding number M (a, b). 

The number M(a,b) depends on the initial data. Thus, (see Corollary on page 45 [3]) under hypothesis of 
the theorem, for any finite value of the integral of 1/ R(t) on [to, oo) there is an open set U of initial data in 
Hi x Hi such that if (Dip,i[>) € U, then there exists a global wave map taking the Cauchy data (ip,ip). In 
particular, this is true for the de Sitter model of universe with R(t) = exp(Ai), A > 0. 

In Robertson- Walker geometry of positive curvature the space time interval is (see Sec. 9. 8 [15], and 



Under assumption that cosmological constant is exactly zero the function a(t) can be found via parameter 
r\ as follows, 



which implies for the function a = a(t) a periodic dependence on time. Although it is not known whether 
the universe actually has exactly periodic behavior suggested by this solution of the Einstein equations, we 
find important to investigate influence of the periodic behavior of metric tensor on the nonlinear waves and 
wave maps propagating in this kind of spacetime. 

Among publications on the time periodic solutions of the Einstein's field equations we shall mention very 
few of them. The study of the periodic solutions to the Einstein's field equations was initiated by Papapetrou 
in [ini [TT1 [TS]. Einstein's old question concerning the existence of solutions free of singularities, vanishing 
and becoming Euclidean at infinity is generalized in |17) for the case of periodic solutions free of singularities. 
The result is negative; such solutions do not exist. It is also shown by Gibbons and Stewart in [6] that any 
asymptotically flat spacetime which is empty and periodic in time is stationary in a neighborhood of null 
infinity. 

Dafermos 0] proved a theorem about the non-existence of spherically symmetric black-hole spacetimes 
with time-periodicity outside the event horizon, other than Schwarzschild in the vacuum case and Reissner- 
Nordstrom in the case of electromagnetic fields and matter sources of a particular kind. This result generalizes 
the so-called "no-hair" theorem from the static to the time-periodic case. That paper addresses the issue of 
the existence of periodic solutions in general relativity in a non-analytic setting, in particular, in a setting 
compatible with the evolutionary hypothesis. For the equations of evolution, time-periodic or stationary 
solutions often correspond to the late time behavior of solutions for a large class of initial data. In the 
general theory of relativity, time-periodic black hole solutions, if they exist, seem to provide reasonable 
candidates for the final state of gravitational collapse. (See [H [9] [10] for more references.) 

For the nonlinear scalar waves and the scalar wave map type equations the influence of the periodic 
behavior of metric tensor on the nonlinear waves and wave maps propagating in this kind of universe was 
found in [24] [25]. It was discovered in those papers that, the generated by periodic coefficient parametric 
resonance interacting with the nonlinearity, in general, leads to the blow up of the solution for arbitrary 
small initial data and for any dimension of the space of spatial variables. 

In order to establish similar result for the wave maps we first prove (Lemma I4.3[) that, in the case of 
n = 2 the only stationary wave maps are the constant wave maps. They evidently have vanishing integral 
(|3j) . This is why we are interested in the small perturbations of the constant wave maps (with the periodic 
in time metric g). 

It was proved in [TTJ Lemma 2.1] that for I = 2, n = 2 and Minkowski space with metric the image 
of the spherically symmetric wave map with smooth initial data belongs to a bounded subset of H 2 . In 



p.131 0) 



ds 2 = -dt 2 + a 2 (t) (d X 2 + sin 2 X i.d8 2 + sin 2 Odcj) 2 )) . 



a = a*(l — cost;), t = a*(?7 — sin 77), 



a» = 



2GM 
3tt 
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particular, for the second component of the wave map we have 

|| lnu 2 |j Lt oo L ~ < oo . (6) 

The main result of the present note is the following theorem that requires some technical assumption 
(Assumption ISIN, see Section [5] for details) on the ordinary differential equation generated by the periodic 
metric. Without loss of generality, we can assume that 

i?'(0) = 0. (7) 



Theorem 1.3 Let R = R(t) of |3|) be periodic, positive, non-constant, smooth function satisfying assumption 
ISIN, and let be Q). Denote (Mf,h) the Riemannian manifold, that is, the half-plane Mf :— {(u l ,u 2 ) £ 

M 2 | u 2 > 0} equipped with the metric tensor hABdu A du B = . ((du 1 ) 2 + (du 2 ) 2 ) . Then, the local solution 

(ir)' 

obtained by small, in the sense of the following integral 

f (\ Qiul i 2 r 57 " 2 1 2 \ j 

J{o}xi» 70+T1+12+ ... +7n=1 G \l( u )' J Y\ u ) 1 J / 

10=0, 1 

some smooth compactly supported spherically symmetric perturbations of the initial data of the constant wave 
map (V,g) — > (M.f,h), cannot be extended globally in time if I £ [0,2). 

If I = 2 and n = 2, then for some smooth compactly supported, with the arbitrarily small integral (0), 
spherically symmetric perturbations of the initial data of the constant wave map, there exist positive numbers 
Co and S such that for all m S N either 

inf \nu 2 {x,m) < —c$e Sm or sup \nu 2 (x,m) > coe Sm . 

Thus, in the case of I = 2 and n — 2, due to the periodicity of the metric of the spacetime, the 
problem in the neighborhood of the stationary spherically symmetric wave map is unstable in the following 
sense: for some smooth compactly supported, with the arbitrarily small integral (j8]), spherically symmetric 
perturbations of the initial data of the stationary spherically symmetric wave map, the estimate ([6|) does 
not hold. 

The proof of the theorem is based on the construction, for every given positive integer number G, of the 
wave map, which takes arbitrarily close to the constant initial data but the image u(tb p ,Xb P ) of some point 
Xb p € K n in the finite time t^p > appears outside of the target manifold Mf, that is, u 2 (tb p ,Xb P ) < 0. In 
that sense the case of I = 0, which implies a linear system for two components of u(t,x), is not exeptional 
and the component u 2 has a finite life span (becomes nonpositive in the finite time). 

This note is organized as follows. In Section [2] we describe the families of geodesies of the target space. 
This geodesies we use to construct wave maps via linear wave equation. In Section[3l we analyze the structure 
of the wave map equations in the Minkowski and Robertson- Walker spacetimes. Section |4] is devoted to the 
finding of stationary wave maps. In Section [5] in order to prove an exponential growth of L^-norm of the 
solutions of the scalar wave equation with periodic coefficients, we derive explicit representations of some 
increasing solutions of the ordinary differential equations with periodic coefficients. In the final Section [6] we 
complete the proof of the main result. 

2 The Target Space. Geodesies 

Let Mf := {(u , u 2 ) \ u 2 > 0} be a half-plane equipped with the metric 

hj du i d^ = r 4v (i du1 ) 2 + ( du2 ) 2 ) » 
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where I G R is a real number, while the metric tensor is 

{h A B) A>B = 1>2 = ( Q x j - 

If I = 2 this is the Poincare upper half-plane model. The upper half-plane H endowed with the hyperbolic 
metric with I = 1 is discussed in Ch.3 [3T]. The Christoffcl symbols are 



rr 1 i = __Lf° 1 

\ AB) A,B=1,2 2u 2 \ l 

The Gaussian curvature is 



i /-l 



A,B=1,2 



2u 2 



1 



/ 



2\Z-2 



Hence, in this family of Riemannian manifolds only the Poincare upper half-plane model and the Euclidean 
half-plane have constant curvature —1 and 0, respectively. 
The geodesies (u 1 ,^ 2 ) = (u 1 (s) , u 2 (s)) satisfy equations 



cPu 1 I du 1 du 2 



ds 2 



i 2 ds ds 



0, 



d 2 u 2 I du 1 du 1 



I du 2 du 2 



ds 2 2u 2 ds ds 2u 2 ds ds 



du 1 



0, 



where s is a natural parameter. If we denote V± := — ; — and V 2 '■— u 2 > 0, then 

ds 



dv L _ L dv 2 = 

ds V 2 ds ' ds 2 ' 2V 2 



I fdV, 



2V 2 V ds 



From the first equation we have 



Vi = C (V 2 ) 1 , C = const e 
By plugging this into the unite speed relation 



(V 2 



we obtain 



with the constraint 



(V2) 



(v 2 y~c 2 (v 2 f = ((v 2 ) s ) 



i(V 2 ) s ) 2 , 



C 2 (V 2 (t)) L < I 



Thus, the function V 2 — V 2 (s) > solves the following equation 



^. = ±^(v 2 y-c 2 (v 2 f. 



Consequently, the system 



^- =±^{u 2 ) l ~C 2 (u 2 f 

d -^ = C(u 2 ) 1 
ds 



has the solution 



(9) 



i\t)=u 1 (0), u 2 {t) = 



(u 2 (0))^ ± 



2-1 



if C = 0, Z 7^ 2, 



■5 



or 



and 



u 1 (s)=u 1 (0), ([u 2 (s 



2-1 



if C = 0, 1^2, 



**(*)= «*(()), u 2 (s) = u 2 (0)e s if C = 0, Z = 2, 



that is a vertical open half-line in the positive half-plane. For C one can eliminate the variable s and 
rewrite the system as a single equation 



du 1 



±^{u 2 ) l -C 2 {u 2 f 



Then, one can integrate it 



± 



C(u 2 ) 



1/2 



^l-C 2 (u 2 ) 



■.du 2 — I du 1 . 



If we denote u := u 1 and v := u 2 , then 



± 



Cx 1 ' 2 



Vl - C 2 x< 



Ax = U — C\ 



We make change of dummy variable x = v t with Ix dx = v dt in the integral, 



Cx 1 / 2 , 2+1 C 

: dx — V 



l Jo y/l - (CV)* 

To evaluate the last integral we use formula (10) Sec. 2. 1.3 2 with 



:dt. 



1 2+1 31 

a=-, b=—, c=- + 1 , c-6-l = 0, z = C 2 v 



l T(b)T(c-b) _ 21 



r(c) 



2 + 1 



and obtain 



2 2 + 1 



21! 1 3 1 2 , 

"^r'2 ; 2 + T ;C v 



where F(a, b;c;() is the hypergeometric function (See, e.g., [2].). In fact, 



(w 1 - Ci) 5 



In particular, for / = 1 we obtain, 



1 



2 + 1 



2( nl 2\2+l 



F 



2+1 1 3 1 



21 ' 2' 2 I 



(u 1 - Ci) 2 = \ -C^/u 2 -C 2 (u 2 ) 2 + arcsin (cV^ 2 ) 



< u 2 < C" 



For / = 2 the geodesies are vertical lines and the upper half-circles given by equation 



(u'-dY + iu 2 ) 2 



1 



d = d 



1 



C 2 ' ~ ' C 
The upper half-circles can be written via the parameter s as follows 

1 e 2s - 1 9 , % 2 e s 



u\s) =d + 



u 2 (s) 



C e 2s + 1 ' 

Thus, we have obtained the second family of geodesies 



C e 2s + 1 ' 



s € 



(10) 



(11) 



G 



To find out the function u 2 — u 2 (s) = v(s) we use the first equation of the system @ and, similar to the 
derivation of (TTU1) . obtain 



2-1 
2 

2-1 



^cv{^ + _l_ c ^ F ( 1|1 + I ; | + I ; ^ 

/l-CW {»'-i(0) + ^cV^OJJ (l, 1 + i; | + i; CV(0)) 



which for the positive u 2 (0) with C 2 (w 2 (0)) < 1 defines an implicit function u 2 = u 2 (s). Then one can use 
(|10[) to find out the function u 1 — u 1 (s). 

The composition of a geodesic with any real-valued solution ip = if)(x, t) of the free wave equation 



1.91 



generates the wave map u = u~ t (see 20J). Such special solutions have been used in [T] to prove, for instance, 
that solution map is not locally Lipschitz continuous. Let 7(s) be an arbitrary curve with values in H 2 , 
which is written in the local coordinates by 7(5) = (71 (s), 72(5)) an d let v(x,t) be an arbitrary real-valued 
function, then the function u(x,t) = 7(^(0;, t)) solves the wave map equation, as soon as v(x,t) solves the 
wave equation and j(s) is a geodesic curve. Thus, we have proved the following statement. 

Theorem 2.1 Let function tp = ip(x,t) be a solution of the covariant wave equation in (V, g) 

g^ (9>-r^)=0. 
Then the following pairs (u 1 ^ 2 ) of functions, 



u 1 (x,t)=Ci, u 2 {x,t) — 



(C 2 )^ ±^<p(x,t) 



, ieR, lei", 2^2, 



u 1 (x,t)=C 1 , u 2 (x,t) = e v ^*\ 1 = 2, 

I e 2(p(x,t) _ 2 2 e f(%,t) 

uHx,t)=d+- e2v{xt) + i , ^ ( ^ t) = ____ ) i eM;XeR « ; l = 2i 
and the function (u (x,t),u (x,t)) — {u l (x, t), v(x, t)) defined by 

u\x,t) = Cl ± ^Cv^F (H+i, I; (| + 1) ; CMM)') , 
^Vl-^V + ^C^F (l, 1 + i; § + I; CV) } 

V^l-CV(O) {^-5(0) + -^ 2 «^ (0)F (1, 1 + j; I + y; CV(0)) } = ^z, t) 



2-; 

are t/ie wave maps u : V — > H 2 . 

3 The Wave Maps in the Robertson- Walker Spacetimes 

In the metric 

g = -dt 2 + R 2 (t)<7, 
in the local coordinates the wave map equation ([T]) reads 

-d 2 u A - nR^R'dtU* + R~ 2 A 7 u A - Y A c d t u B d t u c + R^o^Y^d^ d 3 u c = , 
where we denote by 

A 7 = - a^d k 



7 



the Laplace operator on the manifold with metric <j v} . The wave map solves the system 



I 



-&V - nR^R'dtu 1 + i?" 2 A 7 u 1 + — dtt^dtu 2 - R- 2 a ij ^d^d-v? = 



d 2 t u 2 - nR- l R'd t u A + R~ 2 A 7 u 



;(d t u x dtv} - d t u 2 dtu 2 ) 
-R- 2 a i3 -^{diU l djU l - drfdju 2 ) = 0. 



2w 



If = <5 lJ then we obtain 



dfu 1 + nR-'R'dtu 1 - R~ 2 Au L - —d t u l d t u 2 + iT 2 — VX ■ V x u 2 = 



d?u 2 + nR^R'dtu 2 - R- 2 Au 2 + -J—ldt^dtu 1 - d t u 2 d t u 2 ) 

2u z 



(12) 



-R- 



'^(iv^r-iv^n-o. 



From now on we say that the wave map (it , u 2 ) has a finite life span (blows up in the finite time) if for some 
point (x , t ) with t > the image (m 1 (xo, to), u (xq, t Q )) appears outside of the target manifold H 2 , that is 
either u 2 (xo,to) = or u 2 (xo,to) — oo. 

To reveal the blowup mechanism we use the first wave map of Theorem 12.11 that is, we set in the last 
system u 1 = C\ = const. This corresponds to the choice of the geodesic 7(s), which is the open vertical 
half-line. The second equation of the system reads 



dfu 2 + nR^R'dtu 2 - R~ 2 Au 2 



2u 2 



{{d t u 2 ) 2 ~R- 2 \V x u 2 \ 2 )=i). 



Consider the case of Z 7^ 2. Then we can rewrite the last equation as follows: 



d 2 u 2 + nR-'R'dtu 2 - R- 2 Au 2 - M 1 

1 _ fi-1 

2 ~~ u 



Here 



((d t u 2 ) 2 -R- 2 (\V x u 2 \ 2 )=0. 



7^0. 



2-1 



In particular, if I = 1, then u = 2, while for I = 4 we obtain a = —1. If we denote u :— u 1 and introduce a 
new unknown function v by means of equation 



a 



with parameter a£l and the constants fi, 1/2 = (/x— !)//•*, then from the equations (I12[) we obtain 



u tt + nR^R'm - R~ 2 Au -la[—v + j3) [u t v t - R~ 2 V x u ■ V x v] = 0, 



v tt + nR^R'vt - R~ 2 Av + — f -v + 

2a \fi 



1-2/i 



[(u t f-R- 2 \V x u\ 2 } =0. 



For the case of I ^ 2 and u 1 = const we obtain that the function v solves the following linear equation 

d 2 v + nR- l R'd t v - R~ 2 Av = . (13) 

If we prove that at some point (xo, to) with to > the function v{x, t) takes value —fi/3/a, that is, v(xo, to) = 
— /Lt/3/a, then at that point the wave map (li^u 2 ) blows up. Indeed, it is evident from the definition of v 
that if fi < then at that point u 2 (xo,to) — 00. If fj, > 0, then at that point we obtain u 2 (xo,to) — 0, and 
the image of the point (xo,to) is outside of the manifold Mf. The last case is regarded as blow up as well. 
Note that the half-lines u 1 — const, u 2 > 0, are geodesic in the target manifold H 2 . 
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In the case of the Poincare upper half-plane model H 2 we have I — 2 and 

1 0\ 1 . / 1 AB\ ,..2m( \ 11 



rrM 1 (° M (r 2 1 i /-i o 

\ L AB) A . B = i. 2 - ~2 1 n ' \ as J 



A,B=1,2 ^2^1 0/' ^^;a,B=1,2 u 2 V 1 

and the corresponding wave map equation is 

aV + nR^R'dtu 1 - R^Au 1 - -^(«W)(%* 2 ) + iT^VzU 1 • \7 x u 2 = 0, 
d 2 u 2 + nR^R'dtu 2 - R- 2 Au 2 + \{{d t v}) 2 - (d t u 2 ) 2 ) - Rr 2 ^^^ 2 - \V x u 2 \ 2 ) = 0. 
If we set 

u := u 1 , u 2 = e v > , 

then we arrive at the system 

u tt + nR^R'ut - R~ 2 A 7 u - 2u t v t + 2R~ 2 V x u ■ V x v = 0, 
v tt + nR~ l R'v t - R~ 2 A 7 v + e~ 2v [{u t f - iT 2 |V x u| 2 ] = 0. 



(14) 



If we consider case with u = const, then the second equation implies that v solves linear equation (HU1) 
which has a global solution. Thus, if the initial data for u 1 are the constants, then a global solution of the 
non-linear equation (wave map) exists. On the other hand, if u 2 = const then the first equation is linear, 
and, consequently, it has a global solution. These two families of lines form semi-geodesic parametrization 
of the target manifold. 

4 Stationary Solutions 

We look for the stationary solutions of the wave map system 

7 7 

dfu 1 + nR^R'dtu 1 - R^A^u 1 - — dtv}d t u 2 + iT 2 — V x u • V x u 2 = 0, 

(15) 

d 2 u 2 + nR^R'dtu 2 - R- 2 A y u 2 + -^—(d^dtu 1 - d t u 2 d t u 2 ) ~ R- 2 ^{\V x u x \ 2 - \V x u 2 \ 2 ) = 

2u z 2u z 

with the positive function u 2 (t, x) > 0. The stationary wave map solves the following system of quasilinear 
elliptic equations 

AX - -^V.u 1 • V x u 2 = 0, 

U (16) 
A 1 u 2 + ^{\V x u 1 \ 2 ~\V x u 2 \ 2 ) = Q. 

The second unknown function is assumed to be positive that allows to invoke the Liouville theorem for the 
superharmonic functions. The following statement is evident. 

Lemma 4.1 For every I £ [0, 2] the set of stationary solutions for wave map equation is independent of the 
choice of the function R — R(t). 

Thus for the case of I = 2 and n = 2 we can appeal to the next lemma, which is due to [IT] . 
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Lemma 4.2 f|llj) The image of the wave map belongs to a bounded subset o/H 2 . More precisely, we have 

,1 



\nu 2 \\ LrL ~ < oo, H^IUfZ- <oo. 



(17) 



The bounds depend (at most) on the size of the support as well as some norm ||h[0]||#-i+s, S > 0. 



Lemma 4.3 Suppose that n = 2 and < I < 2. The only stationary solutions of the system HI 5)) with the 
finite integral 



E 



daU 1 



d a u 2 



(u 2 y/2 



dx 



(18) 



and with the positive u 2 (t,x) > are the constant solutions, u (t,x) = c\, u 2 (t,x) = C2 > 0. If 1 = 1, then 
the only stationary spherically symmetric wave maps with target H 2 and finite integral H18\) are the constant 
wave maps. 



Proof. The stationary solution of (IT5|) . that is solution independent of t, v}(t,x) = it 1 (a:), u 2 (t,x) = u 2 (x), 
solves the system (fT6| . In the case of I = we have two harmonic in M 2 functions with the finite integral (fT8|) . 
The Liouville theorem (see, e.g. [El Theorem II]) for the positive function u 2 (x) implies u 2 (x) = const > 0. 
Then, the harmonic functions d a u 1 {x), a = 1, 2, belong to L 2 (R 2 ) only if ■u 1 (a;) = const. 
For the case of I € (0, 2) we set 



u := u , 



2 m 



where = 2/(2 — /), and, consequently, the system (fT6|) reads 



Au 



V x m • V x v = 0, 



AH^Au + - 1)^- 2 |V^| 2 + — (IV^I 2 - ^ 2 v 2 ^ 2 \V x v\ 2 ) = 



Then we obtain 



2(/i-l) 
Au V x u • V x w = 0, 



Av 



+ ^-\V x u\ 2 v 1 - 2 » = 0. 



For the case of < ^ < 2 we have 1 > 0, and, consequently, the second equation of the last system implies 

Av = -^-\V x u\ 2 v 1 - 2 >* <0. 

Thus, the function w = u(x) is superharmonic. According to the properties of superharmonic functions (see, 
e.g., |19j). for the positive solution v = v(x) the last equation implies v = const > 0. Then, the first equation 
of the system implies u = u 1 = u l (x) is harmonic in R™ function. The harmonic function has a finite integral 
(|18[) . that is, the derivatives of this function belong to L 2 (M. n ), only if u 1 (x) = const. 

For I = 2 we set u := u 1 , u 2 = e v and then use (fT4")) . Hence we obtain the stationary solution of the 
system 




Due to Lemma 14.11 and Lemma 
v(x) — lnc > solves equation 



we have v(x) > lnc, where c > 0. The non- negative function w(x) 



AjW - 



^2w 



1 1 2 



0. 
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Then we repeat the arguments have been used in the proof of the case of < I < 2. Lemma is proved. □ 



For n > 2 there is a non-constant bounded superharmonic in 1™ function. (See, e.g., [7].) Moreover, the 
lemma implies the following statement. 

Corollary 4.4 Suppose that < I < 2. If the stationary wave map (V, g) — > (Hj , h) with the finite integral 
H18\) is non-constant, then n > 2. 



This lemma motivates our interest to the perturbations, which have the finite integrals of type (|18|). of 
the constant solutions. 



5 Parametric Resonance in ODE 

Next we make the partial Liouville transformation which eliminates the first derivative Vt in (|13j) with 
(3(x) = const. More precisely, we set 



v = R 2 w . 



then 



Vtt - R~ 2 (t)Av + nR'R^vt = R'% 



w tt - R-\t)Aw + [~R"R- 1 + f (i - \ ) {R' fR- 2 } w 



Thus, we have to study the following linear hyperbolic equation 



w tt - R 



l {t)Aw + { 



--tf'R- 1 



n / n 
2V-2 



(R'fR- 2 } 



w = 



with the 1-periodic positive smooth function R — R(t). 

To this end we are going to apply the Floquet-Lyapunov theory for the ordinary differential equation 
with the periodic coefficients. Consider the ordinary differential equation: 



Wu + {\R- 2 (t) lR"(t)R-\t) - \ g - l) (R'(t)) 2 R- 2 (t)} W 







with the periodic positive smooth non-constant function R — R(t) and parameter A £ 1. 

It is more convenient to rewrite this equation by means of the new positive periodic function 



a(t) =R~ 2 {t), R(t) = {a(t)y 1/2 , 



then 



\a(t) 



3 / a'(t) 
2 V a(t) 



a"(t) 
a(t) 



a(t) 



W = 0. 



Consider now the equation 

ytt(t) + (\a{t)-q{t))y{t) = 
with the periodic coefficients a(t) = R~ 2 (t) and 

"3 (a'(t)\ 2 a"{t) 



(19) 



a(t) 



a(t) 



n tn 
8 V2 



a'(t) 
a(t) 



The first part of the last expression is the so-called Schwarz derivative for the antiderivative of a(t). For 
equation (fT9"l) the spectrum of the eigenvalue problem 



y(0) = = 



11 



is discrete. 

The equation (|19[) can be written also as a system of differential equations for the vector-valued function 
x(t) = t (w t ,w): 

|x(t) = A(t)x{t) , where A{t) := ( \ Q Xa{t) + q{t) \ . 
Let the matrix-valued function X\(t,to), depending on A, be a solution of the Cauchy problem 

±X = A(t)X, X(i ,*o)=( J J ) . (20) 

Thus, W\(i,io) gives a fundamental solution to the equation (|19p . In what follows we often omit subindcx 
A of X\(t,to). The Liouville formula 

W(t) = W(t ) exp QT ^(r)dr^) , 



where := detA(i,i ), S'(i) := Ylk=i Akk[t) w ith ^(*) = guarantees the existence of the inverse 

matrix X\(i,io) • For the matrix X(1,0) we will use a notation 

This matrix is called a monodromy matrix and its eigenvalues are called multipliers of system (|20[) . Thus, the 
monodromy matrix is the value at i = 1 (the "end" of the period) of the fundamental matrix X(t, 0) defined 
by the initial condition X(0, 0) = I (i.e. the matrizant), and the multipliers are the roots of the equation 

det [X(l, 0) - fil] = . (21) 

Due to Theorem 2.3.1 [5] there exist the open instability intervals. We make the following 
Assumption ISIN: There exists the nonempty open instability interval A C (0, oo) for equation 119\) . 
One can find in [5J [H] detailed description of functions a — a(t) and q = q(t) satisfying this condition. 
For instance, in Theorem 4.4.1 [5] one can find asymptotic formula, which allows to estimate the length of 
the instability intervals of the equation obtained from (|19p by Liouville transformation. Then, according to 
next lemma one can find in the instability interval A a number A such that a non-diagonal element of the 
monodromy matrix does not vanish. Moreover, this property is stable under small perturbations of A. 

Lemma 5.1 ([23 -[25 ) Let R(t) be defined on M non-constant, positive, smooth function which is 1- 
periodic. Then there exists an open subset A C A such that &21 7^ for all A G A . 



Next we use the periodicity of b = b(t) and the eigenvalues Mo > 1, Mo" 1 < 1 of the matrix X\(1,0) 
to construct solutions of (Tl9|) with prescribed values on a discrete set of time. The eigenvalues of matrix 
X\(1,0) are ^ and ^q 1 with 6n + 622 = Mo + Mo • Hence {b\\ — Mo) + (022 — Mo) = -Mo + Mo implies 
I &11 - Mol + I&22 - Mo| > I (fen - Mo) + (622 - Mo)| = Imo - Mo *l > °- Tnis leads to 

max{|6n - Mo|, l°22 - Mo|} > tjImo - Mo" 1 ! > °- 
Without loss of generality we can suppose 

l&il - Mol > ^iMo - Mo" 1 ! > °> l°22 - Mo" 1 ! > ^iMo - Mo *l > °> 
because of 6n — Mo = — (022 — McT )• Further, 



1 — =1 — r j— = (Mo - Mo h — r • 

Mo - °22 Mo - 011 o 22 - Mo 
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Lemma 5.2 ([25 ) Let W — Wit), V — V(t) be solutions of the equation 

w tt + (Xa(t) - q{t))w = 

with the parameter A such that 621 7^ and &22 7^ Mo ■ Suppose then that W — W(t) takes the initial data 

W{0) = 0, W t (0) = l, 

and that V = V(t) takes the initial data 

F(0) = 1, V t (0) = 0. 
Then for every positive integer number M G N one has 

W(M) = -JZL^Qtf-^"), 
Mo - Mo 

V(M) = -Mf^4^+Mo- M 7 .^ ia - y 

(Mo - Mo ) (Mo - ohXmo - Mo ) 

In order to prove parametric resonance phenomena in L q spaces with q > 2, we need the following 

Theorem 5.3 Let R(t) be defined on K non-constant, positive, smooth function which is 1-periodic. Then 
for every function ip G C^°(M. n ), which is different from the identical zero, there are positive numbers C v , 5 V 
such that for q G [2, 00] the solutions w(x, t) and v(x, t) of the equation 

w a - R~ 2 (t)Aw + {—H'R- 1 + f (l - f ) {R'fR- 2 } w = (22) 
with the initial conditions 

w(x, 0) = <p(x), w t (x,0) = 0, 
v(x,0) = 0, v t (x,0) = (p(x) 

satisfy 

||t£>(aj,m)|U 9(B ») > C v e s - m Vm G N, 
\\v(x,m)\\ Lq(Rn) >C 9 e^ m VmGN. 

The constant S v depends on the diameter of supp tp only. 

In fact, no = Ho{X) , 6y = %(A), i,j = 1, 2, depend on A as well. 

In fact, the functions W and V depend on A as well, W = W(t, A), V = V(t, A). 
Proof. Let [a, b] C A and tp G C 3O (M). Then by Paley- Wiener theorem 

a b 



^ on I := 



n n 

and for the Fourier transforms of the solutions we have 



where |£| 2 G[a,6] for all (el, 



w(x,t) := (2tt)"/ 2 / e^W(t,\m)dt, 
v(x,t) := (2tt)"/ 2 / e"«y(t,A)^)^- 



Consequently, with some C > we have (analogously to the arguments used in |22| ) 
|w(a;,m)| 2 cfo = \w(£,m)\ 2 d£ > \w(£,m)\ 2 d£ 



1 



> / \W(m,\£\ 2 m)\ 2 dt;> min \W(m,X)\ / |^(0| 2 ^ 
7/ *eMl J 1 
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and similarly for v{x, m). Now, if we take into account the cone of dependence, then for q > 2 we obtain 

\\w(x,m)\\ Lq(Rn) > C v {l+m)~^^w(x,m)\\ L 2 {K n ) > C v> e s "' m Vm e N, 
Mx,m)\\ Lq{Rn) >C^l + m)-^r n \\v(x, m )\\ L 2 (Rn) > C v e s - m Vm e N. 

Theorem is proved. □ 
Corollary 5.4 For a// radfe have 

max |t«(x,m)| > C v e 5,fim , max |t>(a;,m)| > C v e s ' l ' m . 

Theorem 5.5 Le£ i?(i) be defined on R non-constant, positive, smooth function which is l-periodic and 
i?'(0) = 0. Then for every function ip £ C^°(W l ), which is different from the identical zero, there are positive 
numbers Cm, 5 V such that for q € [2,oo] the solutions w(x,t) and v(x,t) of the equation 

d'?v + nR^R'dtv - R~ 2 A^v = (23) 



with the initial conditions 



satisfy 



w(x, 0) = <p(x), w t {x, 0) = 0, 
v(x, 0) = 0, v t (x, 0) = ip(x) 



\\w(x,m)\\ LHRn) >C v e s - m VmeN, 
\\v(x,m)\\ Lq(Rn) > C v e 5 * m Vm e N. 

The constant S v depends on the diameter of supp ip only. In particular, for all m £ N we have 

max \w(x,m)\ > C^e^" 1 , max\v(x,m)\ > Cme 8 *™. 



6 The Completion of the Proof of the Main Theorem 



First we consider the case of < I < 2. We construct the wave map with the finite life-span via a solution of 
the wave equation and the geodesic, which is the vertical line. Let u 2 = const > be a second component of 
the constant (stationary) wave map (u 1 , u 2 ). Consider the perturbation of its second initial data u 2 (x, 0) = 0, 
and look for the global wave map (u 1 (x,t),u 2 (x,t)), such that 

% Q {x)+(3\ , poiVi eC °°(R n ), 
fi = 2/(2 - I). We set <po(x) = and (3 = [u 2 ] > 0. The integral (JSJ of the perturbed initial data 



{0}xR 2 



E 



70=0. 1 







2 

+ 


' d~<u 2 ' 






_(«2)i/a_ 




_(«2)J/2_ 





dx 



E 



{0}xR2 T1+T2 + ... + 

70 = 1 



d^u 2 



T„<G-1 vv ' 



E 

71+72 + - .. + 7n<G-l 



(gVi(g))a(g)^ 

(u 2 y/ 2 



dx 



can be chosen arbitrarily small by appropriate choice of ip± € C^°(]R™). 

Then, by the uniqueness, v}{x,t) = u , and for any real a ^ the function v(x,i) defined by 



v(x, i) = - ( [u 2 (x, t)] 1/M - /?) , u 2 (x, t) > 0, 
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is a solution of the linear equation ([212]) and takes initial values 

v(x, 0) = 0, vt(x,0) = <px(x) ■ 
Hence, for alH > and all x £ W l we have 

u 2 (x,t) = ( -v(a;,i) + (3 

According to Theorem 15 . 5 1 with q — oo and initial condition (|23p. for the appropriate choice of a there exists 
a point (tbp,Xf,p) £ such that 

—v(t b p,x b p) + P < 0. 
M 

Thus, the solution blows up in finite time. 
In the case of I — 2 and n — 2 we choose 

u 2 (x,0)=u 2 , u 2 {x,0)=u 2 ip 1 (x), tp^C^iW 1 ). 

These new initial data have small integral (J5J provided that the function <p\ £ C^°(M. n ) has sufficiently small 
L 2 (M™)-norm. Next we solve equation (|23l) with the initial conditions 

v(x,0)=\nu 2 , Vt(x, 0) = <pi(x) . 

Due to the uniqueness, the vector valued function (u 1 , e vi -- x ^) is a wave map. Theorem 15.51 with q = oo, for 
the appropriately chosen function ipi completes the proof of the main theorem. □ 

Remark 6.1 The perturbation function ip\ £ C^°(R™) can be chosen spherically symmetric, arbitrarily small 
in every given space J- D C^°(R"). 

Acknowledgments 

This paper was originated in November of 2009 when the second author was visiting Osaka University. 
He expresses his gratitude to the JSPS Grant-in- Aid for Scientific Research (A) for the financial support of 
the visit. 

References 

[1] P. D'Ancona, V. Georgiev, Wave maps and ill-posedness of their Cauchy problem, in: New trends in 
the theory of hyperbolic equations, 1-111, Oper. Theory Adv. Appl., 159, Birkhauser, Basel, 2005. 

[2] H. Bateman, A. Erdelyi, Higher Transcendental Functions, v. 1,2, McGraw-Hill, New York, 1953. 

[3] Y. Choquet-Bruhat, Global wave maps on Robertson- Walker spacetimes, Modern group analysis. Non- 
linear Dynam. 22 (2000), no. 1, 39-47. 

[4] M. Dafermos, On "time-periodic" black-hole solutions to certain spherically symmetric Einstein-matter 
systems, Comm. Math. Phys. 238 (2003), no. 3, 411-427. 

[5] M.S. P. Eastham, The spectral theory of periodic differential equations, Scottish Academic Press, Edin- 
burgh and London, 1973 

[6] G W. Gibbons, J.M. Stewart, Absence of asymptotically flat solutions of Einstein's equations which are 
periodic and empty near ininfinity, in: Classical General Relativity. Cambridge: Cambridge University 
Press, 1984, 77-94 

[7] D. Gilbarg, N.S. Trudinger, Elliptic partial differential equations of second order. Reprint of the 1998 
edition. Classics in Mathematics. Springer- Verlag, Berlin, 2001 




15 



S. W. Hawking, G.F.R. Ellis, The large scale structure of space-time, Cambridge Monographs on Math- 
ematical Physics, No. 1. London-New York: Cambridge University Press, 1973. 

D. Kong, K. Liu, Time-periodic solutions of the Einstein's held equations I: general framework. Sci. 
China Math. 53 (2010) no. 5, 1213-1230 

D. Kong, K. Liu, M. Shen, Time-periodic solutions of the Einstein's held equations II: geometric sin- 
gularities. Sci. China Math. 53 (2010), no. 6, 1507-1520: Time-periodic solutions of the Einstein's field 
equations III: physical singularities Sci. China Math. January 54 (2011) no. 1, 23-33 



J. Krieger, Stability of Spherically Symmetric Wave Maps. arXiv:math/0503048vl 



J. Krieger, Global regularity of wave maps from R 2+1 to H 2 . Small energy, Comm. Math. Phys. 250 
(2004) no. 3, 507-580. 

J. Krieger, W. Schlag, Concentration compactness for critical wave maps, arXiv:0908.2474vl 

W. Magnus, S. Winkler, Hill's equation. Interscience Publishers, New York/London/Sydney, 1966. 

H. Ohanian, R. Ruffini, Gravitation and spacetime. Norton, New York, 1994. 

A. Papapetrou, Uber periodische nichtsinguhe Losungen in der allgemeinen Relativitatstheorie. Ann. 
Physik 20 (1957) 399-411. 

A. Papapetrou, H. Treder, Zur Frage der Existenz von singularittsfreien Lsungen der allgemein- 
relativistischen Feldgleichungen, die Teilchenmodelle darstellen knnten, Ann. Physik (7) 3 (1959) 360- 
372. 

A. Papapetrou, Non-existence of periodically varying non-singular gravitational fields, in: Les theories 
relativistes de la gravitation (Royaumont, 1959), .Editions du Centre National de la Recherche Scien- 
tifique, Paris, 1962, pp. 193-198 

J. Scrrin, H. Zou, Cauchy-Liouville and universal boundedness theorems for quasilinear elliptic equations 
and inequalities. Acta Math, 189 (2002) no. 1, 79-142. 

J. Shatah, M. Struwe, Geometric wave equations. Courant Lecture Notes in Mathematics, 2. New York 
University, Courant Institute of Mathematical Sciences, New York; American Mathematical Society, 
Providence, RI, 1998. 

M. Struwe, Variational methods. Applications to nonlinear partial differential equations and Hamiltonian 
systems, Fourth edition. Ergebnisse der Mathematik und ihrer Grenzgebiete. 3. Folge. A Series of Modern 
Surveys in Mathematics [Results in Mathematics and Related Areas. 3rd Series. A Series of Modern 
Surveys in Mathematics], 34. Springer- Ver lag, Berlin, 2008. 

H. Ueda. A remark on parametric resonance for wave equations with a time periodic coefficient, Preprint, 
2011. 

K. Yagdjian, Parametric resonance and nonexistence of global solution to nonlinear hyperbolic equations, 
Proceedings Workshop "Partial Differential Equations", Villa Gualino, Torino, 8-10 Maggio 2000, Edited 
by L. Rodino, Dipartimento di Matematica Dell' Universita di Torino, (2000) 157-170. 

K. Yagdjian, Parametric resonance and nonexistence of the global solution to nonlinear wave equations. 
J. Math. Anal. Appl. 260 (2001) no. 1, 251-268. 

K. Yagdjian, Global existence in the Cauchy problem for nonlinear wave equations with variable speed 
of propagation, in: New trends in the theory of hyperbolic equations, 301-385, Oper. Theory Adv. 
Appl., 159, Birkhauser, Basel, 2005. 



16 



